LOAN COPY: RETURN TO 
AFWL (WLIL-2) 
kirtland AFB, N MEX 


AN ANALYSIS OF 
TERRAIN BIAS ERROR IN 
PLANETARY RADAR ALTIMETERS 


by Richard F. Harrington and William D. Stanley 

Langley Research Center 
Langley Station , Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • MAY 1969 


TECH LIBRARY KAFB, NM 



□i^nao 


III 


AN ANALYSIS OF TERRAIN BIAS ERROR IN 

PLANETARY RADAR ALTIMETERS 

By Richard F. Harrington and William D. Stanley 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



I 


AN ANALYSIS OF TERRAIN BIAS ERROR IN 
PLANETARY RADAR ALTIMETERS 

* 

By Richard F. Harrington and William D. Stanley 
Langley Research Center 

SUMMARY 

Terrain bias error in a pulse-modulated radar altimeter with an omnidirectional 
antenna and a leading-edge range tracker has been analyzed for planetary applications. 
This error is a result of the pulse distortion arising from the surface characteristics of 
a planet. The distortion is due to increasingly longer path lengths at larger incidence 
angles and the variation of the radar cross section with incidence angle. 

The backscatter function of D. O. Muhleman was used as the radar cross-section 
model. An impulse response was first derived as a basis for determining arbitrary pulse 
returns by linear system techniques. An ideal square pulse was assumed for the trans- 
mitted pulse, and a computer program was developed for generating pulse returns. An 
ideal bandwidth- limited receiver was then simulated on the computer and used to filter 
the pulse returns. The position of the received pulse leading edge was defined as the 
time of the 50-percent-amplitude point of the pulse. 

In order to assess the measurable error resulting from the pulse returns, an ideal 
square pulse was used as an input to the receiver simulation program. The position of 
the leading edge of the ideal square pulse after filtering was compared with that calcu- 
lated from the pulse return. The difference between the respective leading edges repre- 
sents the terrain bias error due to pulse distortion. Computations were performed for 
various filter types, transmitted pulsewidths, and bandwidths with input data chosen to 
include the possible range of planetary surface conditions. 

Results of this investigation indicated that terrain bias error varies only slightly 
with filter type and number of poles, is approximately constant for bandwidths greater 
than one-half the reciprocal of the pulsewidth but increases significantly as the bandwidth 
is decreased, and is directly proportional to pulsewidth. 


^Associate Professor of Engineering, Old Dominion College, Norfolk, Virginia 
(Consultant at NASA Langley Research Center). 
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INTRODUCTION 


Altitude measurements are required in many aerospace systems. Altitude informa- 
tion is employed in guidance and control systems, in the correlation of science measure- 
ments, and in the determination of the vehicle trajectory. One method of measuring the 
altitude, when the direction of the surface directly below the vehicle is unknown, is to 
employ a pulse -modulated radar altimeter with a broad beamwidth antenna. One draw- 
back of this approach is that energy is reflected from other directions and is received 
by the broad beamwidth antenna. However, the first return or leading edge of the received 
pulse represents the shortest distance to the surface or altitude. Variations in altitude 
and in surface reflection characteristics distort the received pulse and produce an error 
in the position of the leading edge known as the terrain bias error (ref. 1). 

Two major effects contribute to the pulse distortion. The first effect is associated 
with a nonzero antenna beamwidth in conjunction with the spherical surface geometry. 
Energy leaving the transmitter with increasing look angle from the normal will experience 
increasingly longer time delay and, thereby, a spreading effect is produced. This effect 
can be calculated precisely if minor variations in the surface smoothness are neglected. 

The second effect is the variation of the average or mean radar cross section with 
increasing incidence angle from the surface normal. The variation in radar cross sec- 
tion is difficult to determine precisely for a given planetary body because of the statisti- 
cal nature of the random scatterers and the difficulty of obtaining measurements from 
remote planets. Several mathematical models have been proposed to describe the radar 
cross section of the planetary surfaces (refs. 2 to 5). Although many of these models have 
been found to yield satisfactory results, the Muhleman model (ref. 2) has been used to 
determine the radar cross-section model for both the Surveyor and Apollo lunar module 
programs. Radar cross-section measurements during lunar landing of the Surveyor 
showed good agreement with calculations made with the Muhleman model (ref. 6). There- 
fore, the Muhleman model was used in this investigation to compute the received pulse 
distortion due to variation in radar cross section with increasing incidence angle. 

The purpose of this investigation was to analyze the terrain bias error in a pulse- 
modulated radar altimeter with an omnidirectional antenna. The radar cross-section 
model of Muhleman (ref. 2) was used as a basis for determining the radar returns due to 
an ideal transmitted pulse. The distorted received pulses were computed for both varia- 
tions in altitude and surface reflection characteristics. The effect of filtering the dis- 
torted received pulses in a bandwidth- limited receiver was determined along with the 
effect of filtering an ideal undistorted received pulse. The position of the leading edge 
of the received pulse was defined as the time at which the leading edge had obtained an 
amplitude equal to 50 percent of its maximum amplitude. The terrain bias error, 
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determined by computing the additional time delay of the distorted pulses over the undis- 
torted pulse, was obtained for different filter types, filter bandwidths, and transmitted 
pulsewidths . 


SYMBOLS 

For simplicity, the letters f and g have not been modified when expressed as 
functions of time. For example, the quantities f(t) and g(t) are obtained by replacing 
the angles 9 and <p, respectively, by their appropriate functions of time. The same 
convention applies when the time scale is normalized for f, g, P, S, and y. 

A area, feet2 (meters2) 

a ratio of radar altitude to radius of planet 

B bandwidth of bandpass filter, hertz 

c speed of light, 9.836 x 10^ feet/second (2.998 x 10^ meters/second) 

i(9) normalized Muhleman backscatter function expressed in terms of incidence 

angle 

f(t) function obtained from f(0) when 9 is related to time 

f(v) function obtained from f(t) by normalizing the time scale 

G absolute antenna gain 

G 0 maximum antenna gain 

g($) normalized antenna gain function 

g(t) function obtained from g($) when <fi is related to time 

g(v) function obtained from g(t) by normalizing the time scale 

H radar altitude, feet (meters) 

K amplitude factor for pulse echo 
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k order of low-pass equivalent filter 

n index for computer samples 

Pf. transmitted power, watts 

P r received power, watts 

P(t) received signal resulting from transmitted impulse response 

P(v) received signal obtained from P(t) by normalizing the time scale 

R radius of planet, feet (meters) 

r distance between radar and arbitrary point on surface, feet (meters) 

S(t) received signal resulting from transmitted step function 

T two-way time delay between radar and closest point on surface, seconds 

t time, seconds 

tp transmitted pulsewidth, seconds 

t s time between computer samples, seconds 

v normalized time variable 

x distance in excess of H between radar and planet along arbitrary ray, 

feet (meters) 

x(t) transmitted time function 

y(t) received time function before filtering 

z z -transform variable used in computer simulation 

a Muhleman coefficient 

6(t) impulse function 
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X wavelength, centimeters 

a radar cross section, feet 2 (meters 2 ) 

a Q radar backs catter function 

a Q value of ct q at normal incidence 

y planetary angle, degrees 

c p look angle, degrees 

6 incidence angle, degrees 

t dummy time variable 

DERIVATION OF TERRAIN BIAS ERROR 


For the purpose of this investigation, a mathematical model of the radar cross sec- 
tion as a function of incidence angle was required. Earth based radar measurements of 
the Moon and planets have been used to determine models of the radar cross section. 
However, due to the close proximity of a radar altimeter to a planetary surface, the 
geometry used in analyzing its behavior must be different from the geometry employed 
in developing the model from Earth based measurements. The geometry of the altimeter 
problem is shown in figure 1. Models for the radar cross section have been developed 
as functions of the incidence angle 6. The initial portion of the pulse return is due to 
reflections from the area immediately below the spacecraft. In this region, 6 is very 
small. Other portions of the pulse return are due to reflections from surface areas with 
larger incidence angles. The incidence angle may be uniquely related to the delay time 
for a given pulse return. Models of the angular dependence of the backscatter function 
may be derived from these data. This procedure as used with Earth based measurements 
is reviewed by Evans (ref. 4). 


The fundamental radar equation is derived in such texts as 
Barton (ref. 8). The radar equation relates the received power 
power P^ by 

P t G 2 A 2 a 

p _ k 

r ' (4*) 3 r 4 


Skolnik (ref. 7) and 
P r to the transmitted 

( 1 ) 
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where G is the absolute antenna gain, A is the wavelength, a is the radar cross sec- 
tion of the target, and r is the distance between target and radar. The radar cross 
section is related to the actual area A by 

a = a 0 A (2) 

where a Q is the dimensionless backscatter factor (effective area per unit actual area). 

The preceding relationships are adequate for small targets, but they require modi- 
fication for analyzing the exact behavior of a widely dispersed target such as that encoun- 
tered in a radar altimeter. A differential contribution of received power dP r due to a 
differential cross section da may be expressed as 


P t G 2 A 2 da 
dP r = -£ 

(47r) 3 r 4 

The quantity da may be related to a differential section of surface area by 


(3) 


da = a Q dA (4) 

The backscatter function a Q is, in general, a function of the incidence angle. 

Muhleman (ref. 2) derived a model of the radar backscatter function for a planetary 
surface based on geometric optics and statistical theory. The relationship developed by 
Muhleman contains a single parameter a, which is a function of both the planet roughness 
and the radar frequency. The model developed by Muhleman is expressed as 

cos 0 

a 0 = 5 w (5) 

(sin 9 + a cos 6r 


where a Q is the value of a Q at normal incidence ( 9 = 0°). The quantity a is the 
Muhleman coefficient, which represents the effective mean slope of the planet surface. 
The exact variation of a with frequency is not known. However, it has been determined 
(ref. 2) that for the lunar surface a increases with frequency. The value of a is 
lower for a smoother, more specular surface and is greater for a rougher, more diffuse 
surface. Typical values obtained by Muhleman (ref. 2) are a = 0.265 for the Moon and 
a = 0.1 for Venus, both values being obtained at X = 12.5 cm. 


Equation (5) may be written in the form 

a 0 = 0^(0) (6) 

where f(0) may be defined as the normalized backscatter function. (Note that f(0°) = 1.) 
Insertion of equations (4), (5), and (6) into equation (3) results in 


P t G 2 A 2 a o f(0) dA 
dP r = 

(47r) 3 r 4 


(7) 
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A very convenient means for characterizing a solution of this equation is by using 
the impulse response. The general form of the received power function P(t) due to an 
impulse excitation is derived in appendix A. The function obtained is 


P(t) = f 


K S 2 (t)f(t) 


NT 


( 8 ) 


where K is a constant defined by equation (A23) and T is the two-way delay time along 
the shortest path as defined by equation (All). The quantities g(t) and f(t) represent, 
respectively, the normalized antenna gain function and the normalized Muhleman back- 
scatter function, both expressed as functions of the two-way time required for a particu- 
lar component of the signal to reach the antenna. The time t = 0 is chosen at the begin- 
ning of the return echo. 


The Muhleman backscatter function expressed in terms of time is developed in 
appendix B. In order to simplify some rather unwieldy expressions, the approximation 
t « T was made. This approximation is quite good for most of the cases analyzed. 
With this assumption, the approximate impulse response is given by 


P(t) 


K 

T 



( 9 ) 


The antenna gain must be related to time and placed in equation (9) before the response 
due to an arbitrary pulse excitation may be determined. 

In general, the return echo due to any arbitrary transmitted pulse can be deter- 
mined by convolving the transmitted pulse x(t) with the impulse response P(t). For 
the purpose of this investigation, an ideal square pulse was assumed for x(t). With the 
ideal square pulse, it is not necessary to perform a complete convolution. Instead, the 
step response may be determined as a basis for representing the pulse response. The 
step response S(t) is related to the impulse response P(t) by 

S(t) = f P(t) dt (10) 

J 0 


It is shown in appendix C that for an omnidirectional antenna and with the approximation 
t « T, the step response is 


S(t) = 



-\ 



1 


1 2t 

21 

2 

1 + fl + glT 



L V 2/ aj 




J 


( 11 ) 
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Let tp represent the pulsewidth of the ideal square pulse and y(t) represent the 
received signal. The received signal may be expressed as 

y(t) = S(t) (0 < t < t p ) 

y(t) = S(t) - S(t - t p ) (t > t p ) (12) 

It should be emphasized that y(t) is not the actual RF received signal, but rather it is 
the video (envelope) representation of the signal when a perfect square law detector is 
employed. At the receiver, it is necessary to filter the received signal y(t) with a 
low-pass filter to improve the signal-to-noise ratio. The details of the computer filter 
simulation employed in this investigation are outlined in appendix D. 

Although a pulse return is already distorted when it arrives at the receiver, the 
filter introduces further distortion and time delay. The total delay of a pulse may be 
defined as the time measured to the 50-percent-amplitude point of the pulse. As a means 
of assessing the time delay resulting from the terrain bias, an ideal square pulse may be 
compared with the distorted pulse return. If both pulses are processed through a filter, 
the difference in time delay of the two signals is a measure of the time delay error due 
to terrain bias. Thus, terrain bias error is defined as the difference in time delay, after 
filtering, between a perfect undistorted pulse and the actual distorted pulse return. 

RESULTS AND DISCUSSION 

The parameters that were employed in the investigation of terrain bias error are 
presented in table 1. 


TABLE 1.- DEFINITION OF PARAMETERS 


Variable 

Symbol 

Unit 

Transmitted pulsewidth 

t P 

p sec 

Surface specularity 

a 


Altitude 

H 

feet 

(meters) 

Bandwidth 

B 

Hertz 


Values 

0 . 1 , 1 . 0 , 10 
0 . 01 , 0 . 1 , 1.0 

10 3 , 5 X 10 3 , 5 X 10 4 , 5 X 10 5 

(310), (1.53 x 10 3 ), (1.53 x 10 4 ), (1.53 x 105) 

0.l/t p , 0.3/t p , 0.5/tp, 5.0/tp 


The rationale used in selection of the parameters was to cover values that might be 
employed in a pulse-modulated radar altimeter designed for planetary entry missions. 
The radius of Mars was used in the geometric relationships, and the altitude range of 
1000 to 500 000 ft (310 to 153 000 m) was chosen to represent the range of interest for 
a Martian entry. However, this analysis can be used for any planet by proper scaling of 
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the altitude values by the ratio of planet radii. The values of a were selected to provide 
a broad range that includes the Moon, Mars, and Venus (refs. 2 and 5). A pulsewidth of 
1 p. sec is typical of this type of altimeter, and it was used for most of this investigation. 
Pulsewidths of 0.1 psec and 10 psec were also employed to study the effect of pulsewidth 
variation on the terrain bias error. 

The optimum post-detection bandwidth for the receiver in a pulse-modulated radar 
altimeter is one-half the reciprocal of the pulsewidth (refs. 1, 7, and 8). The bandwidth 
was varied above and below this value to study the effects on the terrain bias error. 

The filter characteristics considered in varying the bandwidth were those of the 
Butterworth (maximally flat -amplitude response), Chebyshev (3-dB ripple), and maximally 
flat-time-delay filters. These filter characteristics are probably the most common types 
encountered in modern radar systems. The differences between the frequency domain 
behavior of the three filters represent a wide variation of trade-off between magnitude 
and phase characteristics. The bandwidth is defined as the width in Hz between the appro- 
priate 3-dB points. The number of poles for each filter was varied from one to four to 
study the effects of increasing stop-band attenuation. For the one-pole filter, all three 
filter types are identical. Thus, a total of 10 different filter responses were considered. 

The pulse responses were computed for the different altitudes and surface speculari- 
ties given in table 1. A pulsewidth of 1 psec was selected for developing most of the data 
in this report. Only the relative amplitudes of the received pulses are presented; that is, 
the pulse amplitudes are normalized to unity maximum value. The pulse response would 
be typical of a video pulse which has been detected by a perfect square-law detector in an 
infinite bandwidth receiver. The normalized pulse responses are presented for various 
radar altitudes and values of Muhleman coefficient in figures 2 and 3, respectively. For 
a = 0.01, the normalized pulse responses at the two lower altitudes are close replicas of 
the transmitted pulse (fig. 2(a)). For a = 0.1, the responses at the two higher altitudes 
approach ramps during the rising portion of the pulse returns (fig. 2(b)). For a= 1.0, 
most of the responses approximate ramps during the rising portion of the pulse returns 
(fig. 2(c)). By comparing figures 2(a), 2(b), and 2(c), it can be seen that increasing a by 
1 decade corresponds closely to increasing the altitude by 2 decades, in terms of the nor- 
malized pulse shape. 

The data for H = 5000 ft (1530 m), 50 000 ft (15 300 m), and 500 000 ft (153 000 m) 
from figure 2 are cross-plotted in figures 3(a), 3(b), and 3(c), respectively, to show the 
effect of a at several different altitudes. 

The pulse responses shown in figures 2 and 3 assumed an infinite bandwidth receiver. 
Radar altimeters require a filter in the receiver to limit the noise input. However, these 
pulse responses are distorted by the use of a finite bandwidth filter. The results of fil- 
tering a 1-psec pulse which has been reflected from a planet's surface with a specularity 


of a = 0.1 in a three-pole Butterworth filter are shown for different altitudes in fig- 
ure 4(a). The amplitude of the filtered pulse is normalized to the maximum value of the 
nonfilter ed pulse echo. Thus, the loss in pulse amplitude to filtering can be determined 
from figure 4(a). The nonfiltered pulse echo distortion for a = 0.1 and a pulsewidth of 
1 p.sec is shown in figure 2(b). A comparison of figures 2(b) and 4(a) reveals that the 
varying pulse distortion with altitude prior to filtering results in a variation of time delay 
with altitude after filtering. This variation in time delay is one of the sources of the ter- 
rain bias error. 

The relationship between pulse distortion due to variation in surface specularity and 
pulse distortion due to variation in altitude has previously been discussed and illustrated 
in figures 2 and 3. Therefore, a variation in surface specularity results in a similar vari- 
ation of time delay of the filtered pulse response and contributes to the terrain bias error. 

Similar normalized filtered pulse responses for the three-pole maximally flat-time- 
delay and three-pole Chebyshev (3-dB ripple) filters are shown in figures 4(b) and 4(c), 
respectively. 

The terrain bias error was computed for a 1-/J. sec-wide pulse for surface specu- 
larities of a = 0.01, 0.1, and 1.0 by using a three-pole Butterworth filter, a three-pole 
maximally flat-time-delay filter, and a three-pole Chebyshev (3-dB ripple) filter. The 
terrain bias error for a 1-p sec -wide pulse and a three-pole Butterworth filter with a 
bandwidth equal to one-half the reciprocal of the pulsewidth is shown in figure 5(a) as 
a function of altitude. Data for constant specularities of a = 0.01, 0.1, and 1.0 are pre- 
sented. Similarly, the terrain bias errors for a three-pole maximally flat-time-delay 
filter and a three-pole Chebyshev (3-dB ripple) filter are shown in figures 5(b) and 5(c), 
respectively. 

Maximum absolute terrain bias error is largest at the highest altitude and for the 
most diffuse surface. A study was made to determine the effect of various filter charac- 
teristics on the terrain bias error under worst altitude and surface conditions (i.e., alti- 
tude of 500 000 ft (153 000 m) and a diffuse surface corresponding to a= 1.0). The 
Butterworth, Chebyshev (3-dB ripple), and maximally flat -time-delay filters were used, 
and the number of filter poles was varied from one to four. The terrain bias error was 
calculated for each combination of filter type and number of filter poles. The terrain 
bias error averaged over all the filter and pole combinations considered was 252 ft (77 m). 
The maximum error was 269 ft (82 m) and the minimum error was 234 ft (71 m). From 
these results, it can be deduced that the terrain bias error varies only slightly with filter 
type and number of filter poles. Therefore, the choice of the filter type and number of 
filter poles would be made from other considerations such as receiver signal-to- noise 
performance, filter construction, and so forth. 



The effect of variation in filter bandwidth and transmitted pulsewidth on the terrain 
bias error was also determined. The bandwidth of the filters was varied from O.l/tp to 
5/tp. The maximum terrain bias error as a function of altitude is shown in figure 6 for 
a 1-ju sec -wide pulse and a three-pole Butterworth filter. The error is essentially the 
same for bandwidths between 0.5/tp and 5/tp and increases for bandwidths less than 

0.5/tp. 

The maximum terrain bias error for a three-pole Butterworth filter with a band- 
width of 0.5/tp is shown in figure 7 as a function of altitude for different pulsewidths. 
The terrain bias error is directly proportional to pulsewidth. Therefore, for a pulse- 
modulated radar altimeter with an omnidirectional antenna and a leading-edge range 
tracker, the pulsewidth determines terrain bias error for an optimum bandwidth and ulti- 
mately provides a limiting factor in the overall accuracy of the radar altimeter. 

The maximum possible terrain bias error due to both variations in altitude and sur- 
face reflection characteristics is one-half the time delay equivalent to a pulse width. 

This error could be reduced significantly if the variation in altitude and/or the surface 
reflection characteristics were known. 

CONCLUSIONS 

Analysis of the terrain bias error of a pulse-modulated radar altimeter resulted in 
the following conclusions: 

1. Terrain bias error is directly proportional to pulsewidth. 

2. Terrain bias error is approximately constant for receiver post-detection filter 
bandwidths greater than one-half the reciprocal of the pulsewidth and increases signifi- 
cantly as the bandwidth is decreased. 

3. Terrain bias error varies only slightly with filter type and number of filter poles. 

4 . The maximum terrain bias error due to both variations in altitude and surface 
reflection characteristics is approximately one-half the pulsewidth. This error could be 
reduced significantly if the surface reflection characteristics were known. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 17, 1969, 

125-22-02-16-23. 
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APPENDIX A 


DERIVATION OF THE RADAR ECHO IMPULSE RESPONSE 


A two-dimensional projection of the planetary-radar-altimeter model is shown in 
figure 1. The altitude H represents the distance of the spacecraft radar above the 
planetary surface, and x represents the additional one-way distance required for a par- 
ticular component of the signal to travel. The angles <p, B, and y are, respectively, 
the look angle, the incidence angle, and the planetary angle. 


Equation (7) in the main body of the report relates a differential element of power 
to the corresponding differential unit of area. The relationship is 


P t G 2 A 2 &of(e) dA 

aP r = 

(4 ir )3 r 4 

The differential unit of area dA may be related to the planetary angle y 
the following equation: 

dA = 2ttR 2 sin y d y 


(Al) 

(see fig. 1) by 
(A2) 


The distance r may be expressed as 

r = H + x 


(A3) 


Insertion of equations (A2) and (A3) into equation (Al) yields 


or 


dP r 

In general, the quantities x 
y. Let 


. PtG 2 A 2 g 0 R 2 f(6>)sin y d y ^ 

2(4tt) 2 (H + x) 4 

and G can be considered to be functions of either <p 
G = G 0 g (<p) (A5) 


where G 0 is the maximum antenna gain, and g(<p) is the normalized antenna gain func- 
tion. Substitution of equation (A5) into equation (A4), and manipulation of the denominator 
yields 


dP r = 


Go 2 A 2 OqR 2 P t g 2 (0)f(g)sin y d y 


2(4tt) 2 H 4 



(A6) 


Now, let 


G 0 2 a 2 QqR 2 

2(4tt) 2 H 4 


(A7) 
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Equation (A6) now becomes 


dP r = 


KiPtg 2 ((ft)f(e)sin y d y 

(* + 1) 4 


(A8) 


A fundamental means of characterizing the nature of the return echo is by using the 
impulse response. Assume then that P^- is a unit impulse — that is, 

P t = 6(t) (A9) 

To utilize the impulse function, it is necessary to relate the differential dy to time. 
Application of the law of cosines yields the following relationship between x and y 
(fig. 1): 

x 2 + 2xH 


cos y=l- 


2(R 2 + RH) 


(A10) 


Assume now that t = 0 is defined as the beginning of the leading edge of the return echo. 
Let T represent the two-way time delay of travel from the spacecraft to the planet along 
the shortest distance. This value is 

2H 


T = 


(All) 


where c is the speed of light. Let r, a dummy variable representing the two-way time 
delay due to the distance x, be given by 

r = ^ (A12) 

c 

and note that 

<A13) 

Manipulation of equation (A10) with the use of equations (All), (A12), and (A13) yields 


H 2 


cos 7=1- 




2r2 ( i + 1) 


For convenience, let 
Use of equation (A15) in equation (A14) results in 


*-i 


(A14) 


(A 15) 


cos 7 


t 


2(1 + a) 


(A16) 
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Differentiation of both sides of equation (A16) yields 


sin y dy = 


a 2 (l + I-j dT 

T(1 + a) 


(A17) 


Substitution of equation (A17) into equation (A8) results in 

a 2 K 1 P t g2(r)f(r) dr 

dPr = 1 t& (A18) 

T(1 + a)(l + g 

where g(<£>) and f(0) have been redefined in terms of t. A contribution of received 
power at time t is due to power transmitted at time t - r. Since the entire time scale 
has been shifted, the received power P(t) can be written as 


P(t) = 


a2Kq P* 
rCi + a) Jq 


g2(r)f(r)6(t - r) dr 


( l + kj 


T(1 + 

A fundamental definition of the impulse function is 

/-»co 

\ h(r)6(t - r) dr = h(t) 

J 0 

Application of equation (A20) to equation (A19) yields 

Pftl E a2Kl fPWW 

' T(1 + ^ (1 + if 

The constant factor may be modified. Let 

a2Ki 


K = 


1 + a 


(A19) 


(A20) 


(A21) 


(A22) 


Substitution of equations (A7) and (A15) into equation (A22) yields 


K = 


; 2>2fv 
o A °o 


2(4tt) 2 (1 + a)H 2 


(A23) 


Equation (A21) now becomes 


P(t) 


K g 2 (t)f(t) 



(A24) 


where the factor T has been retained separately for reasons that will be clear later. 
The function P(t) is the radar echo impulse response. For reasons of convenience in 
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subsequent work, it is desirable to define a normalized time variable v. Let 

v = L = x 

T H 

Substitution of equation (A25) into equation (A24) results in 


P(v) 


K g 2 (v)f(v) 


T 


(1 + v) 3 


(A25) 


(A26) 
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APPENDIX B 


MUHLEMAN LAW APPROXIMATION 


The basic form of the Muhleman equation is investigated, and certain simplifying 
approximations are made. The use of such approximations results in a simpler compu- 
tational form for studying the pulse response. 

The Muhleman backscatter function is given by 


f(e) = 


G^COS 6 

(sin 6 + a cos 0)^ 


(Bl) 


The angle 6 must now be related to time if equation (Bl) is to be employed in the appli- 
cation of equation (A24) or subsequent results. An application of the law of sines to the 
geometry of figure 1 results in 


sin 0 = 


1 + a 


a ( 1 + i) 


sin y 


(B2) 


Substitution of the normalized time variable of equation (A25) into equation (B2) yields 


• n 1 + a 
sin 9 = sm y 

a(l + v) 

The same substitution applied to equation (A16) yields 

a2(y2 + 2v) 

cos y = 1 

2(1 + a) 

Simultaneous solution of equations (B3) and (B4) for sin 9 results in 


sin 9 = 


v + 1 


/v(v + 2) 


1 + a 


a^v(v + 2) 


(B3) 


(B4) 


(B5) 


Generally, almost all of the transmitted energy from an impulse-type excitation will 
be returned within a fraction of the time T measured from the leading edge of the pulse 
return. This condition implies that the normalized time v will be very small in the 
interval of primary interest; that is, v « 1. Furthermore, in the range of altitude mea- 
surements of interest, a = 1. With these assumptions a careful analysis of equations (B4) 
and (B5) yields the following approximations: 

cos 9 ~ — — (B6) 

1 + v 
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sin 6 ~ 


(i + gft + 


1 + V 

Substitution of equations (B6) and (B7) into equation (Bl) yields for the approximate 
Muhleman backscatter function 

(1 + v) 2 


(B7) 


f(v)« 


"'"MX 1 


*%)—■ 

4/ a 


(B8) 


Equation (B8) may now be substituted into equation (A26) to yield, for the impulse 
response, 

Kg 2 (v) 


P(v) * 


T(1 + v) 


1 + 


( i+ !M)XfT 


(B9) 
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APPENDIX C 


DEVELOPMENT OF THE STEP RESPONSE 


The radar echo impulse response was derived in appendix A and simplified in appen- 
dix B. In this appendix, the step response is developed from the impulse response. The 
step response may be used as a basis for determining the actual return from a trans- 
mitted square pulse as indicated by equation (11). 

The step response S(t) is related to the impulse response P(t) by 



(Cl) 


By substituting equations (A25) and (B9) into (Cl), the step response becomes 


S(v) 


- K f 

J 0 


g 2 (v) dv 


(1 + v>£ ♦ (i + |)(i ♦ gjff 


(C2) 


Assume now that the antenna is omnidirectional - that is, g(v) = 1 - and that the region 
of primary interest is where t « T or v « 1. A study of equation (C2) reveals that 
with these assumptions, the step response may be closely approximated by 


S(v) * 


K 

1 + v 




(C3) 


Although the evaluation of (C3) is not completely routine, the details are omitted 
since this operation may be deduced from standard integration procedures. The result 
is 


S(v) = 




(C4) 
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APPENDIX D 


DIGITAL-COMPUTER FILTER SIMULATION 


The process in which a continuous analog filter can be simulated by discrete numeri- 
cal operations on a digital computer is a problem too involved and inappropriately related 
to the principal problem to warrant detailed consideration in this report. Among the 
sources dealing with the details of digital simulation are reference 9 by Kaiser and work 
performed by William D. Stanley at Goddard Space Flight Center while on NASA-ASEE 
Faculty Fellowship in August 1965. A brief discussion of the simulation is presented in 
this appendix. 

The technique employed in the digital simulation was to use a standard linear trans- 
fer function in the p-plane and to map this function to the domain of sampled data system 
(z-plane) by means of the bilinear transformation. This transformation reads 

p = C 1 ~ (Dl) 

1 + z _1 

where z is the z-transform variable and C is a mapping constant. The p-plane trans- 
fer function was assumed to be of the form 


G(p) = 


1 

1 + ajp + a2p2 + 


+ * 40 * 


(D2) 


Substitution of equation (Dl) into equation (D2) yields an equation that becomes increasingly 
unwieldy as k increases. For an arbitrary k, the z-domain transfer function may be 
written as 


, 1 + c.,z - l + c 9 z~2 + . . . + Ci,z - k 

G(z) = -L A ^ 

b Q +b^z - l +b2Z _ 2 + . . . +bkZ - k 


(D3) 


The various constants in equation (D3) may be related to the constants in equations (Dl) 
and (D2) by expansion. The resulting coefficients of digital filters are tabulated for 
k = 1 through k = 4 as follows: 


k = 1 


bo 

1 + cl^ C 

b l 

1 - ajC 

C 1 

1 
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k = 2 

I)q 1 + ftjC + a 2 C 2 
bj 2 - 2a 2 C 2 
b 2 1 - ajC + a 2 C 2 
ci 2 
c 2 1 


k = 3 

b o 

1 + a^C + a2C^ + a3C$ 

b l 

3 + aiC - a 2 C 2 - 3a 3 C^ 

b 2 

3 - a 4 C - a 2 C 2 + 3a 3 C 2 

b 3 

1 - a 4 C + a 2 C 2 - a 3 C 2 

C 1 

3 

c 2 

3 

c 3 

1 

k = 4 


b Q 1 + aiC + a 2 C 2 + a 3 C 2 + a 4 C 4 
bi 4 + 2a x C - 2a 3 c3 - 4a 4 C 4 
b 2 6 - 2a 2 c2 + 6a 4 C 4 
b 3 4 - 2a 4 C + 2a 3 c3 - 4a 4 C 4 
b 4 1 - a 4 C + a 2 C 2 - a 3 C 2 + a 4 C 4 
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The transfer function of equation (D3) may be readily programed on a digital com- 
puter as a numerical algorithm. Let y(n) represent the sampled input to the filter, and 
let y^(n) represent the sampled output of the filter. The computer realization of equa- 
tion (D3) reads 

y f (n) = j^|y(n) + Cjy(n - 1) + c 2 y(n - 2) + . . . + c k y(n - k) 

- b^n - 1) - c 2 y f (n - 2) - ... - c k y f (n - k)j (D4) 

The three types of filters considered were Butterworth, Chebyshev with 3-dB rip- 
ple, and maximally flat -time delay. Detailed consideration of these filter characteristics 
are provided in such synthesis texts as Van Valkenburg (ref. 10) and Weinberg (ref. 11). 
The order of each filter was varied from k = 1 through k = 4. Actually all three filters 
are identical for k = 1. The basic low-pass p-plane normalized function for each filter 
was adjusted to have a reference cut-off frequency of 1 rad/sec. The cut-off frequency 
was defined as the highest frequency at which the response is 3 dB below the direct- 
current level. The coefficients of the reference low-pass analog filters used are as 
follows : 


k = 1 

a l 

Butterworth 

1 

Chebyshev (3 dB) 

1 

Maximally flat-time delay 

1 


k — 2 aj a 2 


Butterworth 

1.4142136 

1.0000000 



Chebyshev (3 dB) 

0.9109423 

1.4125335 



Maximally flat-time delay 

1.3600000 

0.6165333 



k = 3 

a l 

a 2 

a 3 


Butterworth 

2.0000000 

2.0000000 

1.0000000 


Chebyshev (3 dB) 

3.7045854 

2.3832960 

3.9905138 


Maximally flat-time delay 

1.7500000 

. J 

1.2250000 

0.3572917 


k = 4 

a l 

a 2 

a 3 

a 4 

Butterworth 

2.6131259 

3.4142136 

2.6131259 

1.0000000 

Chebyshev (3 dB) 

2.2869936 

6.6056731 

3.2860053 

5.6501357 

Maximally flat-time delay 

2.1300000 

1.9443857 

0.9203426 

0.1960330 
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The constant C in equation (Dl) can be determined from the required mapping 
correspondence between the p-plane and the s -plane. Let p represent the p-plane 
imaginary axis variable and co represent the radian frequency in the s -plane. Assume 
that the time between samples is t s . It can be shown that for co « l/t s , w and p 
are related by 

M «Tpt s (D5) 


Thus, if p r is some reference analog frequency and co r is the corresponding actual 
design frequency, the constant C is 


C = 


2p r 

w r t s 


(D6) 


Accurate simulation of the digital filter was ensured by using 100 sample points 
during the transmitted pulsewidth tp - that is, tp = 100t s . Let B = t where B 
is the bandwidth of the desired filter. With these substitutions and the assumption that 
p r = 1, the constant C is given by 


C = 


100 
77 Btp 


(D7) 
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(a) a = 0.01. 

Figure 2.- Effect of radar altitude on normalized pulse response for values of Muhleman coefficient. 
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Figure 2.- Continued. 
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(a) Three- pole Butte rworth filter. 

Figure 4.- Normalized filtered pulse responses at various radar altitudes, a = 0.1; t p = 1 Msec; b = 0.!j/tp. 
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(c) Three-pole Chebyshev (3-dB ripple! filter 
Figure 4.- Concluded. 
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